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Bose-Einstein condensation (BEC) in a gas has now been achieved [1]. Alkali atoms 
f — . ( 87 i?6, 23 N a and 7 Li) have been cooled to the point of condensation (temperature of 100 

nK) using laser cooling and trapping, followed by magnetic trapping and evaporative 
cooling. This important experimental result has also renewed the interest on theoretical 
studies of BEC [2] . In this contribution we discuss the statistical mechanics of the trapped 
BEC at zero temperature [3]. In particular, we study the stability of the condensate by 
using a variational method with local and non-local interaction between the particles [4] . 

The quantum Hamiltonian operator of a gas of N identical bosonic particles of mass 
m and without spin in a trap is given by 

a ~ 

2-T N h 2 1 N 

H = E ( - |-V 2 + Vo^)) + \ £ V(t u iv) , (1) 

i=l ij=l 

where Vb(r) is the trapping potential and V(r, r') is the interaction potential [3]. Let 
^(r) be the single particle wave function, where a represents a complete set of quantum 
numbers. The totally symmetric many-particle wave function with quantum numbers 
ax, results 



» • p 

where P is the permutation operator and N ai is the number of particles in the state ctj. 
The system is a Bose-Einstein condensate if all the iV particles are in the same quantum 
state a. Thus, for a = ai, we have N ai = N, N a2 = ... = N aN = and the many-body 
wave function reads 

$ a ( ri , r N ) = * a (ri)...* a (rjv) • (3) 
1 Meccanica Statistica di un Condensato di Bose-Einstein Intrappolato. 

2 Presented to SIMAI98, IV National Conference of the Societa Italiana di Matematica Applicata ed 
Industriale, Giardini Naxos, Messina, June 1998. 
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We observe that the mean energy per particle, which is given by 
H 1 f 

e =< — |$ Q >= — J d 3 r l ...d 3 r N <5>* a {r 1 ,...,r N )H<5> a {r l ,...,r N ) , (4) 
can be written as 

e[* a ] = / d 3 v |l|V^(r)| 2 + \/ (r)|M/ a (r)| 2 + ^^ | d 3 rdV|* a (r)| V(r, r')|* Q (r')| 2 , 

(5) 

where e is a functional of the wave function \I/ Q . By considering a weakly interacting gas, 
such that V(r, r') = B5 3 (r — r'), where B = Aith 2 a s /m is the scattering amplitude, the 
mean energy per particle results 

m= /rf 3 r^|V^(r)| 2 + \/ (r)|M/(r)| 2 + ^|M/(r)| 4 , (6) 

where we have omitted the quantum number a and we have written N instead of N — 1 
because we consider a large number of particles. This is the so-called Gross-Pitaevskii 
functional [5]. We observe that the scattering length a s is supposed to be positive for 
87 Rb and 23 Na, but negative for 7 Li. It means that for 87 Rb and 23 Na the interatomic 
interaction is repulsive while for 7 Li the atom-atom interaction is effectively attractive 

[!]■ 

To find the minimum of the energy functional e[^] we impose that the first order 
variation of the energy functional is zero 

M*] = (7) 

with the constraint 



cfr|^(r)| 2 = 1 . (8) 
In this way we get the following Euler-Lagrange equation 

r f) 2 

~ T^V 2 + V (t) + 5iV|^(r)| 2 tt(r) = n*(r) , (9) 

where \i (Lagrange multiplier) is the chemical potential. This equation (Gross-Pitaevskii 
equation) has the form of a nonlinear stationary Schrodinger equation [5]. 

We analyze the Gross-Pitaevskii functional by using a variational technique. The 
external potential of the trap is well approximated by an harmonic oscillator Vo(r) = 
(mu 2 /2)r 2 , where the frequency ui is proportional to the amplitude of the trapping mag- 
netic field. We choose a Gaussian trial wave function for the BEC with a free parameter 
a, which is the standard deviation of the Gaussian, i.e. the mean radius of the condensate 

»« = ^^«p(^)- < 10 > 

The resulting energy function is 

, , 3 , h 2 . 1 3. ma; 2 . 2 BN 1 
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We study when the first derivative of e(a) is zero and we obtain the formula 

For B > there is only one solution, which is a minimum of the energy functional 
(stable solution), instead for B < there are two solutions: one is a maximum (unstable 
solution) and the other a minimum (metastable solution) of the energy functional. In 
fact, for B < the absolute minimum of the energy function is e = oo at a = 0. It follows 



that for B > 0, when iV = we have a = ^Jh/ (ma;) and a — > oo for iV — > oo. For 5 < 
the mean radius a decreases by increasing the number of bosons N to a minimum radius 
a mm , with a maximum number N max of bosons. For greater values of N the condensate 
becomes unstable. It is not difficult to obtain this critical number of bosons. We put 
B = —B and get 

da B 2m J v 2 ; J ' v ; 

from which we have the minimum radius a mm = ^>~ l ^\Jh/ (mw) and the maximum number 
of bosons 

JV^SP^./A. (14) 
5 5 / 4 B 2m\ mu K J 

Thus we obtain an analytical formula of the maximum number N max of bosons for which 

the condensate, with attractive interaction, is metastable. 

We observe that in general we can not neglect the non-zero range dependence of the 

interaction potential. Let us now consider a more realistic non-local interaction between 

the particles 

V(r, r') = B6 3 (r - r') - yr^^ , (15) 

where the constants A, B and T are chosen to reproduce phenomenologically the experi- 
mental data [4]. By using again the trial Gaussian wave function and studying when the 
first derivative of the energy function e(a) is zero we get the following formula 

where erfc(x) = 1 — erf(x) is the complementary error function. For B > there is 
only one solution, which is a minimum (stable solution), instead for B < there are 
three solutions: one is a maximum (unstable solution) and the other two are minima 
(metastable solutions). Thus, we find that with a non-local interaction for B < there 
is a new metastable branch at high density. 

To summarize, we study the stability of a gas of weakly interacting bosons in a har- 
monic trap with local and non-local interaction. To minimize the energy functional, we 
choose a Gaussian trial wave function, where its standard deviation represents the mean 
radius of the condensate. In the case of repulsive interaction the mean radius of the con- 
densate grows by increasing the number of bosons. In the case of attractive interaction 
the mean radius decreases by increasing the number of bosons until it becomes unstable. 
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Our variational technique can be applied also to systems with different shapes of the 
trapping potential. In the future it will be interesting to investigate the time evolution of 
the condensate with exotic initial conditions. 
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